Taste symmetry breaking with HYP-smeared staggered fermions 
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We study the impact of hypercubic (HYP) smearing on the size of taste breaking for staggered 
fermions, comparing to unimproved and to asqtad-improved staggered fermions. As in previous stud- 
ies, we find a substantial reduction in taste-breaking compared to unimproved staggered fermions 
(by a factor of 4-7 on lattices with spacing a « 0.1 fm). In addition, we observe that discretization 
effects of next-to-leading order in the chiral expansion {0{a^p'^)) are markedly reduced by HYP 
smearing. Compared to asqtad valence fermions, we find that taste-breaking in the pion spectrum 
is reduced by a factor of 2.5-3, down to a level comparable to the expected size of generic C'(a^) 
effects. Our results suggest that, once one reaches a lattice spacing of a ~ 0.09 fm, taste-breaking 
will be small enough after HYP smearing that one can use a modified power counting in which 
0{a?) <^ 0{p^), simplify fitting to phenomenologically interesting quantities. 

PACS numbers: ll.15.Ha, 12.38.Gc, 12.38.Aw 



I. INTRODUCTION 

Taste-breaking is a major concern for practical appli- 
cations of staggered fermions. For standard methods to 
work (including the rooting procedure used for dynam- 
ical fermions) taste-breaking must vanish in the contin- 
uum limit. ^ At non-zero lattice spacing, taste-breaking 
is present, and leads to complications in interpreting and 
fitting the results. It is therefore important to find vari- 
ants of the staggered action for which taste-breaking is 
reduced. In this article we study one particularly simple 
approach, hypercubic (HYP) smearing using numeri- 
cal simulations to determine the size of non-perturbative 
taste-breaking, and comparing with other types of stag- 
gered fermions. 

The HYP-smeared staggered action (which we some- 
times call the "HYP action" for short) is obtained by in- 
serting HYP-smcared links in the unimproved staggered 
action. HYP smearing of links was introduced in Ref. 
as a compact smearing approach, and has turned out to 
be a practical and efficacious in several applications. We 
do not repeat the definition of HYP smearing here, except 
to note that key aspects are that it involves multiple (3 
level) smearing and includes reunitarization at each step. 
For more discussion of these and other aspects of HYP 
smearing, see Refs. 0, 0, H, • 

The HYP action is not an O(a^) improved action in 
the sense of Symanzik Q: while taste-breaking O(a^) 
four-fermion operators are removed at tree-level, taste- 
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conserving tree-level O(a^) contributions to the quark- 
gluon vertices and the quark propagator remain. In 
this respect it is inferior to the widely-used asqtad ac- 
tion [1,13], which is tree- level ©(a^) improved. There is 
considerable numerical evidence, however, that the dom- 
inant O(a^) effects are those due to the taste-breaking 
four-fermion operators, and here the HYP action has an 
advantage. The multiple smearing in the HYP action 
leads to couplings between quarks and gluons with mo- 
menta n/a that are suppressed for a larger range of 
gluon momenta and polarizations than in the asqtad ac- 
tion.^ Because of this one expects taste-breaking to con- 
tinue to be suppressed beyond tree-level with the HYP 
action. Indeed, in a recent calculation, Ref. finds that 
the coefficients of taste-breaking four-fermion operators 
induced at one-loop are an order of magnitude smaller 
for the HYP than for the asqtad action. In other words, 
the HYP action turns out to be approximately one-loop 
0{a^) improved for taste- breaking effects. This prop- 
erty, coupled with simplicity, are the main reasons why 
we think it interesting to study the HYP action.'^ 



^ This point is very clearly explained in Ref. @|. See also the 
analysis in Ref. [5]. 

^ Another alternative would be to use the recently introduced n- 
HYP smearing 8]. Here one reunitarizes but does not project 
back into SU(3). This change allows the action to be used in 
dynamical simulations, at no apparent cost to the reduction in 
UV fluctuations. In particular, we expect the reduction in taste- 
breaking to be similar using n-HYP and our "projected-HYP" 
links. Since, however, our ongoing calculations are being done in 
a mixed-action framework (using the freely available dynamical 
lattices generated with the asqtad action) , we do not need to use 
n-HYP links, and it is preferable to complete our studies with the 
same fermion action so as to allow more systematic comparisons 
at different lattice spacings. 
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It should also be noted that one can have both tree- 
level 0{a^) improvement and approximate one-loop im- 
provement of taste-breaking effects. This is accomplished 
by the HISQ ("highly improved staggered quark") ac- 
tion 6]. This action also includes a tuning of the Naik 
term which removes the leading (am)^ error at tree-level, 
thus reducing the discretization errors for heavy quarks. 
This makes it the method of choice for simulating the 
charm quark with staggered fermions. For light quarks, 
however, which are our interest, the lack of complete 
tree-level 0{a^) improvement may not be so important 
in practice, as long as the dominant taste-breaking ef- 
fects are reduced. In particular, the remaining 0{a'^) 
errors are of the same parametric size as those for non- 
perturbatively 0{a) improved Wilson and maximally- 
twisted fermions, both of which are being used for large 
scale simulations. Of course, the key practical issue is 
the relative numerical size of the 0{a^) errors, and we 
return to this point in the conclusions. 

A further advantage of the HYP action is that it avoids 
another well-known problem with unimproved staggered 
fermions — the large size of the one-loop contributions to 
matching factors. For example, perturbation theory for 
Zm, the mass renormalization factor, is very poorly con- 
vergent with unimproved staggered fermions. This is 
solved by any type of link smearing which reduces the 
tree-level taste-breaking four-fermion operators, and thus 
by both HYP and asqtad actions [§] (and, presumably, by 
the HISQ action too). This reduction has been observed 
also for four-fermion operators with the HYP action [l^l . 

With these considerations as motivation, we have un- 
dertaken a thorough study of how HYP smearing reduces 
taste-breaking in non-perturbative quantities. The per- 
turbative improvements described above are important, 
but, if the HYP action is to be phenomenologically useful, 
one must also show that non-perturbative taste-breaking 
is reduced. To do this we study the splittings in the pion 
multiplct. This measure of taste-breaking has several 
advantages. First, very accurate calculations are possi- 
ble, allowing small effects to be picked out. Second, it 
is known from previous studies that the splittings be- 
tween pions are larger than those between states of other 
spin-parities (e.g. the vector mesons or baryons). And, 
third, pion taste-breaking feeds into all quantities, since 
the masses of pions of all tastes enter the non-analytic 
"chiral logarithmic" terms which occur in generic chiral 
expansions. 

It will be useful to recall the pattern of taste-splittings 
within the pion multiplet. The four tastes of stag- 
gered fermions lead to 16 tastes of flavor non-singlet pi- 
ons. These have the spin-taste structure (75 (X) (,f) with 
e {I,£,5,£.^,£.^i5,£.^Li^ = ^Km'Ci']}- For zero spatial mo- 
mentum (which is the case we study numerically) these 
states fall into 8 irreducible representations (irreps) of the 
lattice timeslice group with tastes {/}, {$5}, {^i}, 
{^4}, {Cis}, {C45}, {^ij}, and {^i4}- That with taste $5 is 
the Goldstone pion corresponding to the axial symmetry 



which is exact when m = 0, and is broken spontaneously. 
The properties of the pion spectrum can be studied us- 
ing staggered chiral perturbation theory [l^, It is 
shown in Ref. [l^ that, at leading order in a joint ex- 
pansion in niq and a? , the pion spectrum respects 
an S0(4) subgroup of the full SU(4) taste symmetry. In 
other words, the taste symmetry breaking happens in 
two steps: at 0{p'^) ~ C(a^) the SU(4) taste symmetry 
is broken down to S0(4), while at 0{a^p^) the S0(4) 
taste symmetry is broken down to the discrete spin-taste 
symmetry SW4 As a consequence of this analysis, 

we expect that, to good approximation, the pions will 
lie in 5 irreps of SO (4) taste symmetry, with tastes {/}, 
{C5}, {Cm}' {Cms}' {C^f}- 

A secondary motivation for this work is that the masses 
of pions of all tastes are needed in our companion calcula- 
tion of Bk and related matrix elements. In particular, in 
staggered chiral perturbation theory loop contributions 
from non-Goldstone pions (those with Cf yk^s) is much 
larger than that from the Goldstone pions [ij] . Thus we 
also use this work to test different choices of sources and 
sinks. 

Previous studies have shown that link smearing leads 
to substantial reduction in the splittings in the pion 
multiplet compared to the unimproved staggered ac- 
tion [3, H, [1^. Other features of the action, and in 
particular whether it is fully 0{a'^) improved at tree- 
level, have little effect on these splittings. For example, 
the splittings are reduced similarly by the asqtad action 
(which contains the Naik term) and the smeared p4 ac- 
tion (which contains the "Knight's move" term [l^) [l^. 
The type of smearing used is, however, important. HYP 
smearing (or the double-smearing used in the HISQ ac- 
tion) has been found to be significantly more effective 
than that used in the asqtad action. For example, Ref. Q 
finds that the pion taste splittings'* are reduced by about 
a factor of two when using the HYP (or HISQ) action 
compared to the asqtad action, on quenched lattices with 
a « 0.1 fm. 

Our aim in this paper is to confirm and extend previous 
studies of HYP smearing. We focus on two aspects. First, 
a study of the impact of smearing on the relative size 
of different contributions in the joint chiral-continuum 
expansion. We are particularly interested in the appro- 
priate power-counting to use for valence HYP smeared 
staggered fermions. And, second, a detailed comparison 
with asqtad fermions, using lattices with light dynam- 
ical sea quarks (which are the lattices we are using to 
calculate phenomenologically important quantities). 

This article is organized as follows. The next two 
sections describe some technicalities of our calculation. 
First, in Sec. [Til we explain how we project onto repre- 



^ The quantitative measures used are the spUttings within the pion 
multiplet, A(a^m^), in the chiral limit. We use this measure 
below. 
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sentations of the timeslice group using "cubic sources" . 
Second, in Sec. lIIIi we describe how we construct sink op- 
erators using two different methods, and give examples 
of the fitting of pion correlators. We then turn to our 
numerical results. Sections |TV] and |V] present, respec- 
tively, the pion spectrum calculated using unimproved 
and HYP-smeared fermions on quenched lattices, by 
comparing which we study the impact of HYP smearing 
on taste-breaking, and the nature of the chiral-continuum 
expansion. Section IVII contains our central results: a 
comparison of asqtad and HYP-smeared valence fermions 
calculated on 2 -|- 1 flavor dynamical lattices. We close 
with some conclusions. An appendix discusses the de- 
composition of hypercubic sink operators into timeslice 
representations . 

Preliminary results from this work were presented in 

Ref. [Ulili^- 



II. CHOICE OF SOURCES 

In order to select a specific pion taste we must choose 
sources and sinks that lie in irreps of the time-slice 
group. Here we discuss how this is done for the sources. 
We consider only sources with vanishing physical three- 
momentum. 

We adopt two different methods, which we label the 
"cubic U(l)" and the "cubic wall" sources. Both are 
fairly standard [Ulll^l but, for completeness, we describe 
them briefly. At the end we discuss their relative merits. 

We first fix the configurations to Coulomb gauge. 
Propagators are obtained, as usual, by solving the Dirac 
equation with source h: 

{D + m)x = h (1) 
^ x{x,a;t;A) = ^G{x,a;y,b)h{y,b;A) (2) 

where G is the point-to-point quark propagator, x, y label 
lattice sites, and a, b are color indices. The cubic U(l) 
source at timeslice t is 

h{y, b; t; A) = 5,,,* ^hn+A ^(^' ^) ' (3) 

n 

where n is a vector labeling 2'^ cubes in the timeslice, 
A£ {(0,0,0), (1,0,0),- ■• ,(1,1,1)} labels points within 
the cubes, and r]{n,b) is a U{1) noise vector normalized 
so that 

Ji"^ Yl c') = 5fi^fi'Sc,c' , (4) 

AT— too TV -"^ — ' 

V 

with N the number of random vectors. In practice we 
use only a single noise vector on each configuration. Note 
that the noise vector is the same irrespective of the posi- 
tion in the cube, but differs from cube to cube, and color 
to color. The cubic wall source differs by having the same 



noise vector on all cubes, varying only from color to color: 
hiy, b; t- A) - <5,,,, 5l^_.^^ ry(6) . (5) 

n 

There are eight sources of each type, depending on the 
choice of A. By combining these appropriately, one can 
project onto irreps of the timeslice group. This works 
as follows (using cubic U(l) sources and t = for illus- 
tration). Suppose that we want to select a state with 
spin-taste 75 ® S,f- We calculate propagators from all 
eight sources, a,nd then combine them using the spin- 
taste matrix [23| 

(7s ® (.f)cm ^ ^Tr(7^7s7D7|,) (6) 

(where, following Ref. 7^5 — IC1IC2IC3, etc.,) in the 
following way: 

^ ^' ^' 0; B)e{x)e{A){js ^ ^f) b,a 

■n A,B 

^ E G*{x,a-[0,2n + Alb)in*{n,b) 

V n,A,b 

G(y,c;[0,2m + B],rf)7,(m,d) 

e{x)eiA)ijs&)^F)s,A 
= Y G*{x,a;[0,2n + A],b)G{y,c][0,2n + B],b) 

n,A,B,b 

e{x)t{A){-fs®iF)B,A 
= J2 G{y,c- %2n + B],b){ls®^F)B,A 

n,A,B,b 

G{[0,2n + A],b;x,a). (7) 

The factors of e{x) — (—1)^'^ are needed to hermitian- 
conjugate the second propagator in eq. ([7]). As one can 
see, the final result is the insertion of a bilinear op- 
erator at the source, with quark and antiquark fields 
spread out over a 2^ cube, made gauge invariant by the 
choice of Coulomb gauge, and projected onto zero three- 
momentum. The sum over A and B picks out a bilinear 
lying in a particular time-slice representation. 

In fact, this is not quite correct, because the division of 
the timeslice into a particular set of 2^ cubes breaks the 
single-site translation symmetry. A similar issue arises 
for the sinks, and is discussed in more detail in the fol- 
lowing section and in Appendix \^ The upshot is that 
the cubic U(l) source couples most strongly to the pion 
with the desired taste, but has an additional coupling to 
spin-scalar states with different tastes (and which live in 
a different irreps of the timeslice group). 

This issue does not arise for the cubic wall sources [ill, 
[2]| . Here the bilinear that one is inserting at the source 
is constructed by first summing the fields at a particular 
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position in the cubes over the timeshce, and then com- 
bining them into the desired spin-taste. This projects 
onto a true irrep of the timeshce group. 

We recaU that irreps of the time-shce group contain 
two types of states [25[: those which propagate in time 
with an alternating sign (—1)*, and those which do not. 
This is related to the fact that the bilinears onto which we 
project in Eq. ([7|) involve only a sum over a spatial cube, 
rather than the full 2* hypercube. Thus the timeslice 
operator couples to states with spin-taste 7s and to 
their "time-parity partners" with spin-taste 745s 0^45^, 
the latter propagating with an alternating sign. (Here we 
use a shorthand notation exemplified by 745s — 747575.) 
If we are considering the pion channel then 75 = 75 or 
745, and the alternating state is a scalar, created re- 
spectively by 7s = 74 and 1. These "extra" states do 
not present a significant problem — their contribution is 
reduced relative to the desired pion signal because the 
scalars are heavier, and can be separated because of the 
different time dependence. In addition, it can be reduced 
by using appropriate sinks, as discussed below. A final re- 
duction comes from the fact that operators with 75 = 74 
(at zero momentum and with degenerate quarks) do not 
couple to scalar states in the continuum limit since they 
correspond to conserved charges. 

While our single timeslice sources can produce pious 
of all tastes [lH, the strength of coupling to the pion 
is taste-dependent.^ As can be seen from Eq. ([7]) and 
the definition ©, single-timeslice sources imply that 
S'4 = F4. For pious having tastes with i^4 = 1 (i.e. 
tastes ^5, ^4, and ^^4), the corresponding spin ma- 
trix is thus 75, while for those with F4 = (tastes 1, ^i, 
^45 and ^y ), the spin matrix is 7s = 745- We label the 
former tastes LT (for "local in time" ) and the latter NLT 
("non-local in time" — because the corresponding bilinear 
with the preferred spin 75 involves quark and antiquarks 
on different timeslices). In the continuum limit, opera- 
tors with spin 745 have an overlap with the state that 
is suppressed by rfiT^j K compared to that for spin 75, 
where A ^ {qq)/ f'i is a QCD scale that turns out to be 
large. Thus we expect that our sources are less efficient 
at producing NLT pious than the LT pious. 

The two sources we use represent two extremes: the 
cubic U(l) source corresponds to using a local bilinear 
(local at the scale of a hypercube) summed over the time- 
slice, while the cubic wall source gives rise to a bilinear 
composed of quarks and antiquarks which are separately 
spread out across the entire time-slice. The relative ef- 
ficacy of the sources depends on several factors: which 
has a better overlap with the desired pion, and smaller 
overlap with excited pious and unwanted scalars; which is 
more effective at reducing the noise due to the U(l) ran- 
dom vectors; and whether additional noise is introduced 



^ We note in passing that not all irreps can be produced by single- 
timeslice sources llll . but all pion irreps can. 



by the fact that the cubic source is not a true irrep of the 
timeslice group. In the absence of any clear theoretical 
argument favoring one source over the other, we study 
the issue numerically. 



III. CHOICE OF SINKS AND FITTING 

We use two different methods to construct bilinear op- 
erators at the sink: the "hypercube" or HPC method |23l | 
and the "Goltcrman method" In the former, the bi- 
linear operators (at zero spatial momentum) are 

OsAi) =J2J2 xiz+A)ijs®^F)A,BXiz + B) , (8) 



where n labels 2'^ hypercubes straddling timeslices t and 
t-\-l, with their origins at z = (2n, t), while A and B are 
now vectors in the full hypercube. We stress that these 
HPC sinks always run over two timeslices, even when the 
operator is local in time. They are made gauge invari- 
ant by fixing each timeslice to Coulomb gauge, so that 
spatial links are not required, and inserting appropriate 
time-directed links if ^ F4. More precisely, we av- 
erage over the insertion of the time-directed link at the 
position of the quark and of the antiquark. The links 
that are inserted are HYP-smeared when we are using 
HYP-smeared fermions. 

In the Golterman method, we follow the construction 
of Ref. [11]. For example, if ^4 = -^4, so that the bilinear 
is local in time, it is defined as 

Os.Fit) = J2 ^sA^)xix)^s,FX{x) (9) 

X 

n [(1 - -^^,,1) + \S,,-F^\^^]x{x){10) 

*M x(a;) = ^ [x(a; + A) + x(a; - A)] (11) 

where x — (x, i), and Ts^f{x) is a phase factor given 
in Ref. [llj . Gauge invariance is maintained as for the 
HPC operators. Note that these are single timeslice op- 
erators, unlike the corresponding HPC operators, which 
are spread over two timeslices (even though the terms 
contributing do not involve time links). 

If 5*4 ^ F4, the operators necessarily involve two times- 
lices, and are constructed following the prescription of 
Ref. [U'l, which we do not repeat here. Time-directed 
gauge links added as for the HPC operators. 

The HPC method is, perhaps, conceptually simpler, 
but gives bilinears that are only approximate represen- 
tations of the timeslice group, whereas the Golterman 
method has the advantage of giving true irreps [llj. In 
fact, the HPC operators include the Golterman operators 
as the leading term in an expansion in a, with subleading 
contributions containing derivatives and having different 
spin-tastes. This is discussed in more detail in Appendix 
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1X1 One thus expects the HPC sinks to lead to noisier 
results with more contamination from other states, par- 
ticularly for the cubic U(l) sources which have a similar 
coupling to additional representations. We study this is- 
sue numerically in the next section. 

As noted in the previous section, the pions can be di- 
vided into those having "LT tastes" , which can be created 
by sources with the preferred choice 75 =75, and those 
with "NLT tastes", for which we must use sources with 
is — 745- At the sink end, we always choose 75 = 75, so 
that the bilinears used for the LT pions are local in time, 
while those used for the NLT pions are non-local in time. 



We now briefly describe our fitting method and give 
examples of the quality of our data. For the pions with 
LT tastes, the sink and source have the same spin-taste, 
so that the two-point correlation functions have positive 
time-reflection parity. Thus we fit the data to the follow- 
ing form: 
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f+{t) = Zi[cxp(-mi<) + exp(-rni(it-t))] 

+ Z2(-l)*[exp(-m2i) + exp(-TO2(it-t))I12) 

where Lt = 64 is the lattice size in the time direction. 
The parameters Zi and mi describe the pion contribution 
(spin-taste 75<X)^_f), while Z2 and 1712 describe their time- 
parity partners (with spin-taste 74 C>5 ,^45_f). In this study 
we are concerned only with the pion spectrum, and so do 
not include excited states in the fit function, but rather 
begin fits at times large enough that the form is 
adequate. The typical fitting range is 10 < t < 20, with 
the source a.tt — 0. We use uncorrelated fits (i.e. keeping 
only the diagonal components of the correlation matrix) , 
with errors estimated by the jackknife method. 

We show an example of the data for LT pions in Fig. [1] 
The effective mass mi (t) is obtained by equating the form 
([T^ to the lattice data for times t — {t + 3). The result 
of a fit to the range 7 < i < 16 is also shown. A good 
plateau is seen, with errors of « 1.4%. Note that this 
is a "three-link" pion (quark and antiquark separated by 
three spatial links), which typically has the poorest signal 
of the LT tastes. 

For pions with NLT tastes, our sources and sinks have 
different spins, and consequently the correlation func- 
tions have negative time-reflection parity. The appro- 
priate fitting function is thus 

f-{t) = Zi[exp(-mit)-exp(-mi(it-t))] 

+ Z2(-l)*[cxp(-m2t)-exp(-m2(it-t))] . 



FIG. 1: Effective pion mass for the three- link LT pion, 
am,r(75 ® ^4), for HYP-smeared staggered fermions, cubic 
U(l) sources (at t = 0), and HPC sinks. The fit is dis- 
cussed in the text. Based on 370 quenched gauge configu- 
rations (/3 = 6), and using mi = 7712 — 0.03 (corresponding 
to a particle with mass similar to that of the physical kaon). 
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FIG. 2: Effective mass for four-link NLT pion, am^ (75 (g) 1). 
Parameters are the same as in Fig. [1] 



An example of the corresponding effective mass is shown 
in Fig. [5J There is a reasonable plateau, but with errors 
of « 3% that are larger than those in Fig. [TJ We note 
that this is the pion whose bilinear is the most non-local 
(involving three spatial links as well as the time link) and 
thus has the poorest signal. 



IV. QUENCHED PION SPECTRUM WITH 
UNIMPROVED STAGGERED FERMIONS 

We begin with our results for unimproved staggered 
fermions on quenched lattices. Although, as described 
in the introduction, unimproved staggered fermions are 
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parameter value 

gauge action Wilson plaquette action 

# of sea quarks (quenched QCD) 

P 6.0 

1/a 1.95 GeV (set by p meson mass) 

geometry 16^ x 64 

# of confs 370 

gauge fixing Coulomb 



bare quark mass (am,) 0.005, 0.01, 0.015, 0.02, 0.025, 0.03 

Zm ~ 2.5 



TABLE 1: Simulation parameters for quenched study of unim- 
proved staggered fermions. 

not useful for most practical applications, we begin with 
them because they provide the benchmark against which 
improvements can be measured. Furthermore, we can 
test our methods of calculating the pion spectrum by 
comparing some of our results to those in the literature. 
A final motivation is that we are able to add new infor- 
mation by calculating with a wider range of quark masses 
than used previously. 

The simulation parameters are given in Table [D Here 
Zm is the factor by which the bare quark masses must 
be multiplied to obtain the renormalized mass in the MS 
scheme at a scale of 2 GeV. We have only calculated 
the spectrum for equal quark and anti-quark masses, and 
only consider flavor non-singlet pions (so that there are 
no quark-disconnected contractions). To set the scale, 
we note that am — 0.025 corresponds approximately to 
the physical kaon mass. 

For am = 0.01 and 0.02, we can compare our results to 
those of Ref. [1^ , where the complete pion spectrum was 
calculated using similar methods.^ The comparison is 
shown in Table [TTl Our results are systematically higher, 
although for most pions the difference is only 1 — 2a. For 
the Goldstone, taste <^5 pion, however, the difference is 
larger, ~ 3 — 4(7. These differences could be due to the de- 
tails of fitting (our fitting range starts at somewhat later 
times), and they could at least in part indicate a finite 
volume effect (our lattices are smaller: 16'^ x 64 versus 
20^^ X 40). In any case, we consider the agreement good 
enough for the subsequent discussion of taste-breaking 
effects. 

As a test of our methods, it is also interesting to com- 
pare the relative size of our errors to those of Ref. [1^ . 
The errors are comparable for the Goldstone pion, while 
ours grow more quickly with the distance between quark 
and antiquark in the bilinears. We both use Coulomb- 
gauge sinks (although time-directed links are treated dif- 
ferently), so the differences between the calculations are 



Our results are also consistent with those of Refs. [2ll . l27j but 
we do not show a detailed comparison since both of these works 
used Landau gauge for which the theoretical interpretation of 
the spectrum is not fully justified. 
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FIG. 3: (am^)^ vs. am, for unimproved quenched staggered 
fermions, with cubic wall sources and Golterman sinks. Only 
LT tastes are shown. Solid lines give the results of linear fits 
to the quark-mass dependence. 
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FIG. 4: As for Fig. |3] except using cubic U(l) sources. 



the number of configurations (370 for us, versus 50), spa- 
tial volume (ours is half the size) , and number and type 
of sources (Ref. [2^ uses one wall-source per color, and 
thus three times as many per configuration as we use). 
We speculate that the latter difference, and in particular 
the lack of noise from the U(l) sources, is responsible for 
the smaller errors of Ref. [2g| . 
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Ref. ^ 






This work 


Operator 


am = 0.01 


am = 0.02 


am = 0.01 


am = 0.02 


(75 ® 6) 


0.23877(98) 


0.33476(82) 


0.2430(11) 


0.3382(9) 


(75®Ci5) 


0.2896(14) 


0.3800(12) 


0.2937(25) 


0.3835(16) 


(75 ® Ci4) 


0.3056(16) 


0.3933(14) 


0.3104(27) 


0.3967(21) 


(''Vc (5?) f^") 


3167f2oi 


4016fl6~) 


3230(48 ~) 


4061(28^ 


(75 ® -^') 


0.3300(46) 


0.4110(26) 


0.3323(57) 


0.4127(40) 


(75 'S'Us) 


0.2922(30) 


0.3826(18) 


0.2978(62) 


0.3866(35) 


(75 ® ^ij) 


0.3085(26) 


0.3957(18) 


0.3108(46) 


0.3980(33) 


(75 ® Ci) 


0.3168(33) 


0.4026(19) 


0.3211(49) 


0.4061(38) 



TABLE II: Comparison of results for (am^) with unimproved staggered fermions with those from Ref. 
those with Golterman sink operators and cubic wall sources. 



Our results are 
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FIG. 5: {amn) vs. amq for NLT tastes using unimproved 
quenched staggered fermions, with cubic wall sources and 
Golterman sinks. Solid lines give the results of linear fits 
to the quark-mass dependence. 



We show the comparison of the two sources for LT 
pions in Figs. [3] and [H and for NLT pions in Figs. [5] and 
ini^ In both cases the resuhs from the two sources are 
consistent. We note, however, that the cubic wall sources 
lead to noticeably smaller errors for NLT pions than the 
cubic U(l) sources, although there is no such difference 
for LT pions. 

We show results with HPC sinks (and cubic wall 



''' We do not focus here on the detailed chiral behavior, and thus use 
linear fits. As the figures show, these give a good representation 
of our data over the range of quark masses. We are not concerned 
that the Goldstone pion mass does extrapolate exactly to zero at 
vanishing quark mass, since this can be explained by the missing 
higher-order terms in the chiral expansion. 



FIG. 6: As for Fig. Oexcept with cubic U(l) sources. 



sources) in Figs. [7| and [51 Comparing these with Figs. [3] 
and [5l respectively, we see that the choice of sinks has 
no significant impact on the signal. (This is shown in 
more detail for HYP fermions below.) Thus we conclude 
that the extra irreps coupled to by the HPC sources do 
not significantly degrade the signal. Nevertheless, all else 
being equal, we prefer the Golterman sources. 

We now turn to the nature of the taste breaking with 
unimproved fermions. Having results for all the tastes 
over a range of relatively light quark masses (down to 
~ jj^phys^g-j a^Q.y^fg dlseutangle the different con- 

tributions. We note that the splittings are compara- 
ble to the pion mass-squareds themselves at the light- 
est quark masses, indicating that the appropriate power- 
counting is 0(0^) « 0{p'^). As seen in previous stud- 
ies, the states are ordered according to the number of 
gauge-links traversed when moving between quark and 
anti-quark fields in the bilinear, although the splittings 
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are not directly proportional to this number. In addi- 
tion, we clearly observe that the slopes differ between 
the tastes: the splittings all increase with mj, with al- 
most a factor of two increase between the chiral limit 
and am — amP^'^^ w 0.025. (The values of the slopes 
for LT tastes are given in Table IVll below.) This differ- 
ence in slopes arises at 0{a^p^) in chiral power-counting, 
and thus is a next-to-leading-order (NLO) effect. We 
thus expect it to be a small correction for our lightest 
quark masses, increasing to a significant correction for 
m ~ m^^^^. This is indeed what we find. 

Staggered chiral perturbation theory predicts that the 
NLO contributions to the slopes for the four different 
LT tastes are independent [23|. In fact, the numerical 
results tell us that the slopes for the heavier three tastes 
are consistent with each other, while being significantly 
larger than that for the Goldstone pion. This implies 
that some of the coefficients in the NLO staggered chiral 
Lagrangian are significantly smaller than others. 

Another NLO prediction is that SO (4) taste symme- 
try should be broken, i.e. that the leading order (LO) 
degeneracy between tastes Ai and A4, Tn and T^j, and 
V4 and Vi, will no longer hold. The splitting is caused 
by terms of ©(a^p^), and thus affects only the slopes 
{arriT^Y /{aruq), and not the intercepts [11]. The differ- 
ence in slopes within SO (4) irreps is predicted to be of 
the same order as that between SO (4) irreps. Taken at 
face value, this would imply that the S0(4) symmetry 
would be completely broken at the point that NLO con- 
tributions approach the size of those of LO (which, as we 
have seen, occurs at the upper end of our mass range). 

In fact, we find no such breakdown in the SO (4) sym- 
metry. The accuracy to which this symmetry holds can 
be seen most easily from Table [TTl and also by compar- 
ing Figs. [7| and [5] (the points with the same colors in the 
two figures live in the same SO (4) multiplet). The er- 
rors for NLT tastes, while larger than for the LT tastes, 
are small enough that we can rule out a scenario with 
complete breakdown of SO (4) symmetry for our largest 
quark masses. In fact, our results are consistent with 
SO (4) symmetry for all quark masses. Thus we conclude 
that our results are consistent with the prediction of no 
SO (4) breaking in the intercepts, and that the coefficients 
of the SO (4) breaking in the slopes must be much smaller 
than the generic size. 



V. QUENCHED PION SPECTRUM WITH 
HYP-SMEARED STAGGERED FERMIONS 

Using the same set of quenched gauge configurations as 
in Sec. IIVI we now study the pion spectrum with HYP- 
smeared staggered fermions. The parameters we use are 
summarized in Table Hill We recall that HYP smearing is 
straightforward to implement in practice: one simply ap- 
plies HYP smearing to the gauge configuration and then 
calculates the propagators using the unimproved stag- 
gered action on the resulting lattice. HYP smearing is 
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FIG. 7: (am^)^ vs. amq for LT tastes using unimproved 
quenched staggered fermions, with cubic wall sources and 
HPC sinks. 
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FIG. 8: As for Fig. [7] except for NLT pions. 

carried out after Coulomb gauge fixing. 

Despite using larger bare quark masses, our calcula- 
tions with HYP-smeared fermions correspond to phys- 
ical quarks which are slightly lighter than those used 
for unimproved staggered fermions in the previous sec- 
tion. This is because Zm ~ I for HYP-smeared staggered 
fermions, whereas k, 2.5 for unimproved staggered 
fermions. In particular, the physical strange quark mass 
is amPj^^^ ~ 0.052, slightly larger than our largest quark 
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parameter 



value 



gauge action 
# of sea quarks 
/3 

1/a 

geometry 
# of confs 
gauge fixing 
smearing metliod 
bare quark mass {amq) 



Wilson plaquette action 
(quenched QCD) 
6.0 

1.95 GeV (set by p meson mass) 
16^ X 64 

370 
Coulomb 
HYP (II) of Ref. [10] 
0.01, 0.02, 0.03, 0.04, 0.05 
fa 1 



TABLE III: Simulation parameters for the quenched study 
with HYP-smeared staggered fermions. 



mass [23], whereas the physical strange quark mass for 
unimproved fermions is approximately equal to the next- 
to-largest quark mass in that calculation {amq — 0.025). 

Table IIVI shows the dependence of pion masses on 
source and sink type for one choice of quark mass. As 
for unimproved quarks, the sink choice is unimportant 
(the central values are consistent, and the errors are very 
similar). It also remains true that the cubic wall sources 
lead to smaller errors for the NLT pions, although the 
reduction in the errors is a much smaller effect. While 
we do not understand why this is, it does mean that, 
in practice, the choice of source appears less crucial for 
HYP-smeared staggered fermions. 

The resulting pion mass-squareds are shown for our 
preferred source and sink, and for all tastes, in Fig. [9l 
Comparing to the previous figures, we see that taste- 
breaking is substantially reduced. * A quantitative mea- 
sure of this is given in Tabic IVl which collects the mass- 
squared splittings after linear extrapolation to the chiral 
limit, both for the HYP and the unimproved action. Note 
that this measure is insensitive to the particular range of 
quark masses chosen (as long as they are small enough 
that linear fits work well). For the LT tastes, where the 
splittings with the HYP action are statistically signifi- 
cant, we find that HYP smearing reduces the splittings 
by a factor of 5 — 10. The ordering of states is unchanged. 
For the NLT tastes, the splittings are reduced by simi- 
lar factors, but the errors are large enough that a more 
quantitative statement cannot be made. For the same 
reason, we cannot make any quantitative comments on 
the size of possible SO (4) sphtting. 

Our results for the reduction in taste-breaking are 
qualitatively consistent with those obtained previ- 
ously 0, H, llHl (although differences in simulation details 
do not allow a quantitative comparison). Having the re- 
sults for a range of quark masses allows us to draw some 
conclusions on the appropriate power-counting. In par- 



The vertical scales in the plots for the unimproved and HYP 
action differ because the quarks, and thus the pions, are hghter 
with the latter. 



ticular, given the smallness of the splittings, it may be 
appropriate for this range of quark masses to treat the 
these O(a^) effects as of NLO rather than of LO, i.e. to 
assume that ^ ^ rn^ . 

A test of this idea is provided by the taste-dependence 
of the slopes in the vs. plots. As noted in the 
previous section, this difference is an 0{a?p^) effect, and 
thus of NLO in the standard power counting of staggered 
chiral perturbation theory. It is observed at the expected 
order of magnitude with unimproved staggered fermions 
(see, e.g., Fig.[3l). If the power-counting changes to ~ 

for HYP-smeared fermions then the difference in the 
slopes would be of NNLO (next to next to leading order), 
and thus expected to be very small for our range of quark 
masses. This is indeed what we find. Table IVII collects 
the slopes for the LT tastes, and we see that, with HYP- 
smeared fermions, they are equal for all tastes within the 
statistical uncertainty, which is at the 2% level. 



VI. COMPARISON OF ASQTAD AND 
HYP-SMEARED STAGGERED FERMIONS 

In this section we compare two types of improved stag- 
gered quark actions: the "asqtad" action used by the 
MILC collaboration and the HYP-smeared action. We 
use one of the so-called "coarse" lattice ensembles gen- 
erated by the MILC collaboration [s^], with the aver- 
age light quark mass roughly one fifth of the physical 
strange quark mass. The parameters of the study are 
summarized in Table IVlTl We note that the conventions 
of the asqtad action are such that, when comparing to 
the HYP action, the appropriate bare quark mass to use 
is aniq/uQ, where uq is the "average link" [3l| . In other 
words, if the HYP and asqtad quark masses are nomi- 
nally equally, the HYP masses are actually smaller by a 
factor of uq — 0.86774 (this value coming from Ref. 31 1 
using one of the conventional definitions of mq). To set 
the scale, the strange quark mass with the asqtad ac- 
tion is « 0.045, i.e. near the upper end of the range 
used, while that for the HYP action is « 0.057, which 
is somewhat above our heaviest quark mass. The results 
for asqtad valence quarks are from Ref. [3l|, while those 
for HYP-smeared quarks are new. 

We only have results for HYP-smeared quarks using 
the cubic U(l) source. Although the results of the previ- 
ous section indicate that this source is slightly inferior for 
NLT pions, our errors are small enough for the purposes 
of this study. 

When we calculate HYP-smeared staggered quark 
propagators using the cubic U(l) source, we choose the 
source time slice randomly configuration by configuration 
such that unwanted autocorrelations in the two-point 
pion correlation functions are removed. 

In Fig.[TOl we show the results from Ref. [sij for (m^)^ 
for all tastes as a function of quark mass for asqtad va- 
lence quarks. This should be compared to our results 
with HYP-smeared valence fermions, given in Fig. II H and 
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HPC 






Golterman 


Operator 


GUI 




cw 


GUI 


GW 


(75 ® 6) 


0.2262(27) 




0.2261(26) 


0.2263(27) 


0.2267(35) 


(75®Ci5) 


0.2313(27) 




0.2311(31) 


0.2314(30) 


0.2302(26) 


(75 ® Ci4) 


0.2355(29) 




0.2348(32) 


0.2356(29) 


0.2339(27) 


(''Vc (5?) fi") 


2395f3oi 




2375f3li 


2396('311 


2374(^29^ 


(75 <8>C45) 


0.2293(74) 




0.2239(69) 


0.2343(72) 


0.2239(69) 


(75 ® Cij) 


0.2338(77) 




0.2264(56) 


0.2346(79) 


0.2233(61) 


(75 «) Ci) 


0.2392(79) 




0.2281(56) 


0.2409(65) 


0.2262(58) 


(75 «> 


0.2407(83) 




0.2298(62) 


0.2403(84) 


0.2309(65) 



TABLE IV: Source and sink dependence of (om,r). Results using HYP-smeared staggered fermions at am — 0.02 on the 
quenched lattices. HPG and Golterman denote the choice of sink, while GUI and GW denote cubic U(l) and cubic wall 
sources, respectively. 
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FIG. 9: (am^)^ vs. aruq for HYP-smeared staggered fermions in quenched QGD, with cubic wall sources and Golterman sinks. 
Linear fits to the quark-mass dependence are shown. The left panel shows LT tastes; the right panel shows NLT tastes. 



Taste (F) 


Af (unimproved) 


Af (HYP) 




0.0233(14) 


0.0025(11) 


Ci4 


0.0320(27) 


0.0050(14) 


^4 


0.0425(49) 


0.0073(13) 


C45 


0.0244(91) 


0.0036(53) 




0.0296(73) 


0.0050(51) 




0.0303(78) 


0.0076(47) 




0.0511(147) 


0.0078(56) 



TABLE V; Gomparison of taste symmetry breaking be- 
tween unimproved and HYP-smeared fermions. Here Ap — 
[am-^{F)]'^ — [arriTriS.s)]'^ , extrapolated to the chiral limit. Both 
results are obtained using Golterman sinks and cubic U(l) 
sources in quenched QGD at /3 = 6.0. 



Fig. [T^l Note that different axes are used from the earlier 
plots, with the scale here being set by ri, which is 2.61a 
on these lattices. 



Taste (F) 


C2 (unimproved) 


C2 (HYP) 




5.57(2) 


2.28(4) 




6.08(4) 


2.26(4) 


Cj4 


6.03(9) 


2.24(4) 


^4 


5.99(10) 


2.24(4) 



TABLE VI: Gomparison of taste symmetry breaking in slopes 
between unimproved and HYP-smeared fermions. C2 is de- 
fined by fitting the data to [am,r(F)]'^ = ci + C2{amq). Both 
results are obtained using Golterman sinks and cubic wall 
sources in quenched QGD at (3 = 6.0. 



Comparing the two valence actions, we see that the 
taste-breaking is substantially reduced for HYP-smeared 
fermions. To make this quantitative we quote the values 
for some of the taste splittings, extrapolated to the chiral 
limit using linear fits, in Table IVini The splittings with 
HYP smearing are smaller by a factor of 2.5 — 3. This is 
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parameter 



value 



gauge action 
sea quarks 
sea quark masses 

P 
a 

geometry 
# of confs 
gauge fixing 
valence quark type 
valence quark mass (asqtad) 
valence quark mass (HYP) 



1-loop tadpole-improved Symanzik 
2+1 flavors of asqtad staggered 
ami = 0.01, arUg — 0.05 
6.76 
0.125 fm 
20^ X 64 
640 (asqtad)/ 671 (HYP) 

Coulomb 
asqtad and HYP staggered 
0.007, 0.01, 0.02, 0.05 
0.005, 0.01, 0.015, 0.05 
?a 1 



TABLE VH: Simulation parameters used for the comparison 
of asqtad and HYP-smeared valence quarks. The HYP smear- 
ing uses the parameters of case (H) in Ref. [ici| . 




2mqri 



FIG. 10; (rim^)^ vs. 2rimq on unquenched configurations, 
with asqtad valence quarks. All tastes are shown. The cor- 
respondence of the taste labels to those in other figures is as 
follows: TTs =PxP, TTos = PxA4, TTis = PxAi, TToi == PxTi4, 
TTij = PxTij ttq = PxV4, TTi = PxVi, and vr = PxS. Results are 
from Ref. I3l|. The lines are linear fits to each taste. 




0.3 



2 niq ri 



FIG. 11: As for Fig. \T0\ for LT tastes except using HYP- 
smeared valence staggered fermions, with Golterman sinks 
and cubic U(l) sources. 




consistent qualitatively with the result noted in the intro- 
duction from Ref. obtained on quenched lattices. It is 
also comparable to the factor of 3.6(5) improvement ob- 
tained using HISQ rather than asqtad fermions on similar 
dynamical lattices [6]. 

A surprising feature of the comparison of Figs. llOl with 
Figs. [11] and [12] is that the slopes for the Goldstone pi- 
ons differ: that for HYP-smeared fermions is about 1.4 
times smaller than that for asqtad fermions. (The dif- 
ference is emphasized by our use of the same vertical 
scales.) This slope is a parameter in the corresponding 



FIG. 12: As for Fig. \W\ for NLT tastes except using HYP- 
smeared valence staggered fermions, with Golterman sinks 
and cubic U(l) sources. 



chiral Lagrangian, and should agree for the two actions 
up to differences in renormalization factors for the quark 
masses (aside from taste-conserving scaling violations). 
In fact, as already noted, the asqtad masses should be 
multiplied by 1/uq « 1.15, so that the horizontal scale 
in Fig. [TU] should be stretched by a factor of 1.15, thus 
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Taste {F) 


Af (asqtaa) 


Af (MYF) 


5*5/^45 


0.205(2) 


0.0708(12)/0.062(13) 




0.327(4) 


0.1380(29)/0.125(16) 




0.439(5) 


0.2008(62)/0.199(26) 


I 


0.537(15) 


0.257(34) 



TABLE VIII: Taste symmetry breaking effect: Af = 
[rimvr(-F)]^ — [f 111-^(55)]^ extrapolated to the chiral limit. 
Note that the scale is set here by ri and not the 1/a used 
earlier. Results from unquenched lattices comparing valence 
asqtad (from Ref. [3l[) with HYP-smeared fermions (Goiter- 
man sinks with cubic U(l) sources). 



reducing the slope by this factor. This leaves a factor 
of about 1.2 difference in slopes to explain. This could 
be entirely due to differences in — the two-loop re- 
sult for asqtad fermions is ~ 1.16 [s^l, while that for 
the HYP action on asqtad sea-quarks is not known yet 
at one-loop. It is important to note, however, that the 
splittings given in Table rvTIIl are in the chiral limit and 
thus are unaffected by any relative factors between the 
bare quark masses. 

Returning to taste-breaking, we note that for both 
types of valence fermion the slopes for different tastes 
are nearly equal. For asqtad fermions we also see that 
no SO (4) breaking is apparent even at the largest quark 
masses (which, we recall, are close to the physical strange 
quark mass). Thus, despite the relatively large splittings 
with asqtad fermions [of 0{a^) in the chiral counting], 
which are clearly to be treated as comparable to the con- 
tributions from the quark mass [of ©(p^)], there are no 
indications of NLO contributions [of 0{a^p^)]. This is 
perhaps surprising and is in any case in contrast to unim- 
proved staggered fermions where 0{a^p^) effects were 
clearly visible for a similar range of quark masses. For 
HYP-smeared fermions there is also no indication of NLO 
effects. The slopes are very close and, as can be seen from 
Table IVHll there is no significant SO (4) breaking. 

We next compare taste-breaking for HYP-smeared 
fermions between the quenched lattices considered ear- 
lier and the dynamical lattices. Clearly, the splittings 
are significantly larger for the latter. To make this quan- 
titative we compare the splittings in physical units in 
Table HXl where we see that they are 2-2.5 larger on the 
unquenched than on the quenched lattices with HYP- 
smeared fermions. Can we understand this? The answer 
is yes, at least qualitatively. Taste-breaking is expected 
to scale as a^a^, since the contributions from one-gluon 
exchange have been removed. The lattice spacing on the 
dynamical lattices is larger, with the ratio of to that 
on the quenched lattices being about 1.5. We also know 
that a at a scale of « 1/a (which is the scale relevant for 
taste-breaking gluon vertices) is larger on the dynami- 
cal lattices, since the long-distance physics has been (ap- 
proximately) matched by construction, and the coupling 
runs toward zero more slowly for the dynamical lattices. 
Combining these effects, we expect to find that the taste- 



breaking for HYP-smeared fermions is perhaps twice as 
large on the dynamical as on the quenched lattices. This 
is in accord with our results. Note, however, that a^a^ is 
reduced by about a factor of 2.5 between the MILC coarse 
and fine lattices |i31i] . Thus we expect the taste-splitting 
on the MILC fine lattices with the HYP-action to be 
similar to that we observed above on the a = 1/10 fm 
quenched lattices. 

In summary, the main difference between asqtad and 
HYP-smeared valence quarks is a substantial reduction 
in the LO taste-breaking splittings between SO (4) mul- 
tiplets. Unlike the quenched case studied above, how- 
ever, the power-counting ~ appears to remain nec- 
essary on these coarse MILC lattices. It seems likely, 
however, that taste-breaking will be a NLO effect with 
HYP-smeared quarks on the fine MILC lattices. 



VII. CONCLUSION 

In this paper, we have compared three versions of 
staggered fermions — unimproved, HYP-smeared, and 
asqtad — from the standpoint of taste-symmetry break- 
ing. We confirm and make quantitative the hierarchy ex- 
pected from previous work, namely that HYP smearing 
is significantly more effective at reducing taste-breaking 
than asqtad improvement. This is summarized in Ta- 
bic IIXI Our results imply that HYP-smeared valence 
fermions on the fine MILC lattices (with a « 1/11 fm) 
should have taste-breaking small enough that it will be 
possible to treat discretization errors as a NLO effect in 
chiral perturbation theory (rather than as a LO effect as 
required with asqtad fermions, and with HYP-smeared 
fermions on the coarse MILC lattices) . This is an encour- 
aging result for our ongoing calculations of electroweak 
matrix elements using HYP-smeared valence fermions. 

To put the results into context, it is useful to con- 
vert the size of the splittings into a scale characterizing 
discretization errors. Taking the splitting of the 2-link 
LT pion (taste ^^j^) as indicative of an "average" split- 
ting, and setting Am^ = a^A"*, we find A = 0.5 and 
0.6 MeV for HYP fermions on the quenched and un- 
quenched lattices, respectively. These are scales within 
the range of sizes that one would expect for generic dis- 
cretization errors. In other words, the improvement due 
to HYP smearing brings the taste-breaking errors down 
into the range expected of generic discretization errors, 
so that discretization errors with staggered fermion are 
no longer anomalously large. 
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Taste(F) 


C[U6Ilch6cl ffi P 

A F (unimproved) 


- 0.1 fm) 

Af(HYP) 


unquenclied (g ^- 
AF(asqtad) 


s 0.125 fm) 

Af(HYP) 


Cms 


0.0931(56) 


0.0094(41) 


0.0759(7) 


0.0262(5) 




0.128(11) 


0.0189(53) 


0.1211(14) 


0.0511(11) 


Cm 


0.170(20) 


0.0279(49) 


0.1625(19) 


0.0743(23) 



TABLE IX: Comparison of taste-breaking for LT tastes on all lattices considered in this paper, 
now in units of GeV^. 
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APPENDIX A: ANATOMY OF THE 
HYPERCUBE OPERATORS 

In Sec.lIIH we defined the hypercube operators, eq. ([8]), 
and noted that they are not irreps of the time-slice group 
of Ref. [ll| . Here, following the procedure of Ref . [33] , we 
express them in terms of true irreps. For simplicity, we 
work in the free theory, i.e. we do not explicitly include 
gauge links, but it is straightforward to add these at the 
end, following, say, the prescription in the text. 

Consider a single hypercube operator of general spin 
and taste: 

Ofp = E xiA)hFWl^)ABXiB) . (Al) 

A,B 

Summing this over a timeslice as in eq. ([5]) leads to the 
HFC operators we actually use, but we delay this sum- 
mation until necessary. To simplify the discussion, we 
first assume that 15"^ — ~ 1 for one index fi, while 
\Si, — — for h' ^ ji. In other words, we first con- 
sider a "one-link" bilinear. We also need the "averaging" 
operator $p given in eq. (jlip and the derivative operator 

D.xiA) - ^[xiA + i,)~xiA-fl)], (A2) 
in terms of which 

x{A + fi) - ^,aiA) + D^x{A) (A3) 
x{A-fi) = <^>,x{A)-D^x{A). (A4) 

Now we can express the single-hypercube operator as 



(summation over A, B now implicit) 



0|f - X{A)ijs^^F)AB{SB,,A, + l+SB,,A,-l)x{B) 

= lx{A){-fs <E) ^F)ABSB,.A,+i['i>,.x{A) + D^x(A)] 
+ li'^f.xiB) ~ D^x{B)]SA,,B,-ihs ®iF)ABX{B) 
+ lx{A)iis ® ^f)abSb,.a,~i['^^^x{A) - D^x{A)] 

+ + D^x{B)]SA,,B, + lhs ® ^F)ABXiB) 

= \x{A){-fs <E) ^f)ab[^^.x{A)] 
+ ^[$pX(S)](7s®ef)ABX(S) 

+ \x{A)h^5S ® ^^5F)AB[Df,x{A)] 

- l[D^x{B)]h^.5S <E> W)asx(S) 

In the last step we have used the following result for spin- 
taste matrices, assuming jS'p — = 1: 



(7S S.f)ab{B - A)t, = (7^55 ® S.t.bF)AB ■ (A5) 

Now we reinstate the sum over spatial hypercubes. 
Then the following identity holds: 

x{A){is®^F)AB[^f.x{A)] 

n 

= X^[$MX(B)](7S®efWx(S) 

n 

= E x{is®^f)x]t (A6) 

n 

where [• • • ]t indicates that this is part of an irrep of the 
timeslice group of Ref. [l3|. In fact, if is a spatial 
index, the [• • -jr part of the HPC operator is identical 
to the corresponding Golterman operator of eq. ([9]). If 
/i = 4, the operator in (jA6p transforms in the same irrep 
as the corresponding Golterman operator, although it is 
not identical. 

To pick out the remaining parts of the HPC operator. 
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we define, as in Ref. [33| 

n 

= J2 {x(^)(7m5S ® W)as[£'a.x(^)] 

n 

-[^mX(5)](7/x5S ® e^5F)Ai3X(S)} • (A7) 

This also transforms as part of an irrep of the timeshce 
group [a different irrep from that in which the operator 
(eq. IA6p transforms, unless /i = 4]. Thus we can express 
the HPC operator as a sum of two true representations: 

n n 

+^IZx(7/.5S<» W)'^mX (A8) 

n 

Note that for = 4 these two operators, while having 
different spin-taste, lie in the same irrep of the timeslice 
group. This is the standard "doubling" of states dis- 
cussed in the text. 

The extension to operators containing more than one 
link is straightforward. For example, for distance two 
operators (i.e. — = \Si, — F^\ = 1 for IJ. ^ v, and 
\Sp - Fp\ — Q for v), we find: 

E^^.F = [Ex(7^^)x]r 



Thus there are four irreps (all different \i p ^ A ^ v) 
contained in the HPC operator. 



,7m5S 






$ ^iv5f) 



In general, the HPC operator includes the Golterman 
operator (or, if 1 15*4 — -^4 1 — 1, an operator in the same ir- 
rep) as the leading term as well as operators with deriva- 
tives having different spin-taste. One expects the lat- 
ter to have a coupling to physical states suppressed by 
^ aAqcD for each derivative. A similar analysis holds 
for the cubic U(l) sources, while the cubic wall sources 
should produce only true irreps aside from the effects of 
using only a few noise vectors. One should also keep 
in mind that the symmetries are not completely imple- 
mented when one has a finite ensemble of gauge configu- 
rations, although we expect this effect to be smaller than 
the use of only one source-set per configuration. Thus 
we expect the contamination from coupling to states of 
unwanted spin-taste to be smallest for the combination 
"Golterman- wall" , followed by "Golterman-U(l)" , then 
"HPC-wall", and finally "HPC-U(l)". This is not, how- 
ever, the only source of systematic problems. The cou- 
pling to excited states within the same irrep should be 
more significant for wall than for U(l) sources. 
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